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Abstract. We study the motion of kinks in a two-lane model of the totally 
asymmetric simple exclusion process with open boundaries. We analytically study 
the motion of the kinks by a decoupling approximation. In terms of the decoupling 
approximation, we find that the positions of the kinks become synchronised, though 
the difference in the number of particles between lanes remains non-zero when the 
rate of lane change is asymmetric. The validity of this result is confirmed for small 
asymmetric cases through the Monte Carlo simulation. 



1. Introduction 

We often encounter congestion of pedestrian flows and traffic flows in our daily life. 
We also observe the stuck of grains in granular flows. It is important to study the 
mechanism of congestion not only from the industrial point of view but also from the 
physical point of view. For the sake of scientific research, we need to analyse a simple 
model which captures the essence of the phenomena. 

The asymmetric simple exclusion process (ASEP) is one of the simple models 
adequate to describe such the transport phenomena [1, 2]. It is a stochastic system 
of particles moving asymmetrically on a lattice. The simplest limit of ASEP is that the 
particle is only allowed to hop in one direction, which is called the totally asymmetric 
simple exclusion process (TASEP). 

It is known that the stationary state of one-lane ASEP under open boundary 
conditions has been obtained exactly [3, 4, 5, 6]. Dynamical properties of TASEP are 
also studied expensively. The exact solutions of the master equation by Bethe ansatz 
on an infinite system [7] and a periodic system [8] have been obtained. Furthermore, 
the current fluctuations in an infinite system and a semi-infinite system [9] are also 
studied. In the open boundary system, we can draw a phase diagram by the parameters 
of inflow rate and outflow rate at the boundaries. On the phase boundary between the 
low density phase and the high density phase, there exists a diffusive domain wall (kink) 
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[10, 11, 12, 13, 14]. Recently, Takesue et al [14] have derived a /"^^^ law in the power 
spectrum as a function of the frequency / based on the random walk picture of the 
kink, and confirmed its quantitative validity from the comparison of the Monte Carlo 
simulation with their theoretical prediction. 

However we little know the properties of a multi-lane ASEP which is more realistic 
than the one-lane model. There are several two-lane models of ASEP [15, 16, 17, 18]. As 
used in [15] or [16], a reahstic lane change rule should refer to the states of the front sites. 
However, the rule makes us difficult to analyse because we need to construct a transfer 
matrix to refer to the states of three or four sites, which are the current site, the front 
site and the side site, or the front site of another lane. Belitsky et al [15] successfully 
analysed the long-time properties of such the two-lane model in an infinite system. Here, 
we adopt a simpler model of lane change in which the particle may change lanes when 
the side site is vacant and do not refer to the front site. Although the two-lane(channel) 
model which adopts this simple rule is dealt by Pronina et al [18] who analysed the 
model based on a cluster approximation and compared the result with their extensive 
simulations, the lane change rates and the boundary parameters are symmetric in their 
model. We extend their model to the case of asymmetric lane change rule and boundary 
parameters to study a more general situations. 

The purpose of this paper is to clarify the motion of kinks in a two-lane TASEP. 
To fulfil the analysis, we introduce our model and explain how to specify the position 
of the kink in the next section. In section 3, we discuss the motion of two kinks based 
on a decoupling (mean field) approximation. We find that the motion of the kinks are 
synchronised though the number of particles in one lane is different from that in another 
lane. We compare the solution with the results of Monte Carlo simulation. In section 
4, we discuss the validity of the mean field approximation. We find that the two-point 
correlation function is small during the relaxation process from the independent motion 
of two kinks to a synchronised motion of them. In section 5, we conclude our results. 

2. Our model 

2.1. Introduction of our two-lane model 

Our two-lane model is defined on a two lane lattice of L x 2 sites, where L is the length of 
one lane. We introduce the occupation variable Tj.^ where Tj.^ = 1 and Tj-i = represent 
the occupied state and the vacant state on the jth site in the £th lane, respectively. 
The particle move forward by the rate 1 during the time interval dt if the front site 
is vacant. We assume that all the particles drift from the left to the right. The open 
boundary condition is characterised by the inflow rate ae and the outflow rate (3^. The 
particle is injected to the system by the rate when the flrst site in the ^th lane is 
empty, while the particle is extracted from the system by the rate Pe when the Lth site 
in the ith lane is occupied. On all sites, the particle is allowed to change lanes only 
to the neighbouring site. A particle on the 1st (2nd) lane can change lanes by the rate 
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Figure 1. A sample picture of the two- lane model. 



rj(r|) when the side site in another lane is empty. 

We denote the probabihty of finding the system in the configuration r = 
{ti;i, T2;i, ' ' ' , 7"l;2} by -P(ti;i, T2;i, • • • , tl-2) ■ We write the time evolution of the two-lane 
model by the master equation, 
J ^ 



2 L-1 



+ XI XI X (^.7J+l;^)-rj;€,r,+i;€;o-3;^,(T,+i;^-P(Tl;l, " ' ' , CTj-i, CTj+l-l, " ' ' , TL;2) 

^=1 i = l + 
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+ X X (^i;l,2)<Tj;i,cTj;2;rj;i,T,;2-P('''l;l, • • • , CTj;!, " " " , Crj;2, • • • , Tl;2) 

where aj-e is used for a dummy variable in the summation, and the transition matrices 
hi-e, hi-e, hjj+i-e, /ij;i,2 are represented as 
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Here the density function (r,-.^) and two-point function {Tj.£Tk-/') are defined by. 



{Tj-ZTk-/') = y^Jj-tTk-l'PjTi-l, ■ ■ • ,TL;2), 



(4) 



where the summation is taken over all the configurations. The time evolution of (r^j^) 
is written as 

d 



dt 



(Tj;^) - Jj-l,j;i - - Jj-/^e' + Jj;i'-^e 



(5) 
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for 7^ where the current Jjj+i-e between site j and j + 1 is 

Jj,j+i;e = {Tj-eC^ - Tj+i-i)) (6) 
and the currents between lanes are 

Jj;l-.2 = n(Tj;l(l-r,;2)) J j -2^1 = r ^ {Tj .2(1 - Tj .i)) . (7) 

2.2. The Position of the Kink 

It is known that a kink appears when the inflow rate is equal to the outflow rate and 
both rates are smaller than 1/2 in one lane ASEP. The kinks also appear in the two-lane 
model when ai = Pi < 1/2 and q;2 = < 1/2. For r| 7^ and 7^ 0, it is obvious 
that the motion of one kink depends on another kink. 

We need to specify the position of the kinks in order to discuss their correlated 
motion. The position of a stable kink in the one-lane ASEP can be determined by using 
the second class particle [10, 11]. However, we adopt another method to determine the 
position of the kink by the whole number of particles in a lane. This deflnition has been 
used in the domain wall theory [12, 13, 14], and give the exact position of the kink when 
the inflow and outflow rates are small. The advantage to adopt this method is that it 
is much simpler than the method by the second class particle. 

We introduce (A^^) for the whole number of particles in each lane 

L 

im = J2i^r,i)- (8) 

We also introduce (No) and (Nr) by (No) = ((A^2) + {Ni))/2 and (A^^) = {N2) - (Ni) 
respectively. The position of the kink xg is defined from the equation based on a kink 
picture; 

xe = M^^, (9) 

Pt,- - Pl;+ 

where pi.^± represent the density of the Ith lane. The index + represents the right side 
of the position of the kink and index — represents the left side of the position of the 
kink. It is straightforward to give the eq.(9) by summing up the equation 

{Tj-/) = Pt- + {Pi;+ - Pt,-)0{j - Xg) (10) 

from = 1 to L, where 6{z) is the step function, 

^ ^ f 1 for z > 
[0 for 2; < 

Thus, once (A^) is known, we can determine the position of the kink. 
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3. Mean-Field Theory 

For the large system size L, we can neglect the boundary terms. Thus the equations for 
(iVc), {Nr) are given by 

^W = (12) 

d 



dt 



{Nr) = - 2t^{Nr) + 2(ri - rt) $^(r,;i(l - r,.,^)). (13) 



When = r^, the eq.(13) is reduced to 

^(Ar«) = -2rt(7V«). (14) 

Thus {Nr) relaxes to exponentially, and the number of particles becomes identical in 
both lanes in the long time limit. 

However for ^ r^, the problem becomes nontrivial because of the two-point 
function in the second term in the right hand side of the eq.(13). In general, the two- 
point correlation function is determined by an equation including three-point correlation 
function. Thus, we cannot obtain the exact form of the many-point correlation functions 
without truncation of the hierarchy of correlation functions. Here we adopt the simplest 
truncation, which is the decoupling (mean-field) approximation as 

Y,{r,,,{l - r,,,)) ^ X^(r,;,)(l - {r,,,)). (15) 

We also use the kink picture (10) to approximate the density profile {rj-^f). 

Let us discuss the motion of two kinks starting from the initial condition where 
two separated kinks exist in both lanes. The density profile changes in time during 
the synchronisation of the kinks. Furthermore, we assume that the density changes 
by keeping the density profile (10). Thus we have to determine the time evolution of 
the density pt,± on the both sides of the kink. By introducing Pg;± = ^'^'^^^^'^ and 
Pr\± = P2;± — Pi;±, the time evolution equations for pG-± and pr-± are respectively 
written as 

l^m. = (16) 

Equations (16) and (17) can be solved exactly, 

PG;± = P°G;± (18) 



where 
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n + ± A/(n + ^t)' - 2(n - ^t)'(«i + «2 - («i + ^2)72) 

"^-^ v;^, ('^^ 

C± = ^^'l^ - (22) 

The initial conditions p^.j^ and p5t;± ^-^s taken as the stationary densities in one-lane 
model as 



PG;- = PK;- = «2 - «i (23) 

= "'T"' = - (24) 



where 



Thus we obtain the density 

7 g~7* ^ e~^^ 

pi;- = pi-7P^F3a: = + 7i=^^^ ^'^^ 

and the density 

7 g~7* ^ e~^^ 

Pi = ^ P2 = ^ (27) 

Therefore we obtain the time evolution of the density profile (t^j^) for rj^ t^. 
Substituting eq.(lO) into (13) with the aid of (15) we obtain, 

j^{Nr) = -{2r^ + e(l - p2;- + Pv,+)){Nr) + 2e{p^,+p2.,-L + (1 - - P2;-)(A^g)) (28) 
for xi < X2 and 

^(Nr) = -(2r| + e(l - P2;+ + Pi;-))(iVi?) + 2t{pi,.p2,+L + (1 - pi;_ - P2;+)(iVG)) (29) 
dt 

for Xi > X2- Equations (28) and (29) can be solved exactly, though the expression is 
lengthy (see (A. 4) and (A. 5)). Here, we present the solution for {Nr) in the long time 
limit as, 

= (30) 

for xi > X2i and 

= !^i(wpi (31, 

- ep'2 

for xi < X2- These results show that there remains the mean difference of the number 
of particles between lanes. The validity of our analysis based on the decouphng 
approximation is confirmed by the comparison of our result with the Monte Carlo 
simulation when |r| — rj^| is not large. Figure 2 shows the quantitative accuracy of 
our analysis in the time evolution of (Nr) . 
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Figure 2. The figure shows the comparison of time evolution between the simulation 
and the calculation. The solution (A. 4) is shown by the solid line and the simulation 
result is shown by the cross(x). The parameters are ai = /3i = 0.1, a2 = P2 = 
0.15, fj = 0.11, r-f = 0.1, L = 1000. Time step is taken for each Monte Carlo step. The 
initial condition is fixed to {Nji)o = 75 and averaged over 1000 samples. 



Though {Nn) remains finite, the positions of the kinks are synchronised. In fact, 
from eqs. (30) or (31) and (9) we obtain, 

X2 - xi = 0. (32) 

Thus the positions of the kinks become identical in the long time limit. This result is 
reasonable, because we cannot choose a preferable congestion front in traffic jams. 



4. Discussion 



Now let us discuss the validity of the decoupling approximation. Although it is difficult 
to evaluate the two-point function exactly, it is possible to evaluate it from the Monte 
Carlo simulation. The result of our simulation for the two-point function 

E^=l Ej = l('^j;^) 

is shown in Fig. 3 for the same parameters as used in Fig. 2, where we realise that A is 
small between t = 30 and t = 100. The synchronisation is realised before i = 30 as we 
can see in the Fig. 2. Thus, we may expect that the decoupling approximation adopted 
here works well to describe the synchronisation of the kinks. 

However there is a certain parameter region that the decoupling approximation 
fails. In fact, the left figure of Fig.4 shows {Nji) obtained from the simulation deviates 
from that in the decoupling approximation. The positions of the kinks are not identical 
in this case. The value A given by the simulation in the right figure of Fig.4 when it 
deviates from the decoupling approximation. The value A is almost in the region 
between t = 30 and t = 200. In such the region the stationary density p£-± is not given 
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Figure 3. The figure shows the time evolution of the value A from the simulation. 
The parameters used in this simulation are the same as those used in Fig. 2. 
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Figure 4. The left figure shows the comparison of time evolution between the 
simulation and the calculation. The solution (A. 4) is shown by the solid line and 
the simulation result is shown by the cross(x). The right figure shows the time 
evolution of the value A from the simulation. In both figures, the parameters are 
ai = /3i = 0.1, a2= 02 = 0.15, = 0.03, = 0.01, L = 1000. Time step is taken for 
each Monte Carlo step. The initial condition is fixed to {Nii)o = 75 and averaged over 
1000 samples. 



by the eq.(17). The comparison of p£-- between the result of the simulation and the 
solution of the eq.(17) is shown in Fig. 5. We can see that the density derived from the 

mean-field theory deviates from the result of the simulation. Thus in the region that 
the final densities of each lane differ from the solutions of the eq.(17), the final positions 
of the kinks are not identical. 

To verify the valid region of the decoupling approximation, in Fig. 6, we plot the 
deviation from the decoupling approximation in the difference between two kinks as 
a function of rj/r-|- — 1. From Fig. 6 we find that the deviation is proportional to 
5^ = (rj/r| — 1)^ for the small asymmetric cases. This is an interesting result, because the 
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Figure 5. The figure sliows tlie difference between the stationary density given by 
the simulation and that by the eq.(17), in the corresponding low density region. The 
plus(+) stands for the first lane, the cross(x) stands for the density in the second 
lane and the approximation is shown by the solid line. The boundary parameters are 
ai = 0.1, /3i = 0.9, a2 = 0.15, (32 = 0.85, = 0.03, r| = 0.01, L = 1000. 
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Figure 6. The figure shows the final difference between the kink positions given by 

the simulation and the analysis. The cross(x) shows the difference and the solid line 
is the guideline which is proportional to (— — 1)^. The boundary parameters are 
ai=pi= 0.1, a2=P2 = 0.15, L = 1000. 



decoupling approximation predicts (Nji) oc 6. Namely, the decoupling approximation 
can be used when we can regard S as finite but 5^ as negligible. In addition, the curvature 
of the parabola is relatively small, which ensures that we may use the approximation 
in 5 < 0.2. The quantitative analysis of the violation of the decoupling approximation 
will be a future problem. 
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5. Conclusion 



In this paper, we have studied the motion of kinks in the two-lane TASEP. We obtain 
the exphcit time evolution function of the average number of particles in each lane which 
is related to the position of the kink by adopting the decoupling approximation of the 
two-point correlation function. We find that the positions of the kinks are synchronised, 
though the number of particles in a lane can be different from that in another lane. 
We confirm the validity of our analysis by comparing the result of the Monte Carlo 
simulation and the analytical result. The deviation from the mean-field analysis is 
small when the lane change rates are nearly symmetric. 
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supported by the Grant-in-Aid for Scientific Research (Grant No. 15540393) of the 
Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan, and 
the Grant-in-Aid for the 21st century COE 'Center for Diversity and Universality in 
Physics' from MEXT, Japan. 



Appendix A. The calculation of (A^ij) 



In this appendix, we give the explicit expression of {Nr). Equation (29) is solved as 
{Nr) = (ArB)oe--^odt'(2n+e(i-P2;++pi;-)) 

+ e rdM*'^*"(2^T+^(i-P2;++.,-))(2Lpi^_P2;+ + 2(1 - p,,_ - P2;+){Ng)) (A.l) 
Jo 

Here, we perform the integral in the argument of the exponential function. 



Therefore, 



r 1 

/ dt' (1 - P2;+ + Pl;-)^ 2 p[t+- In 

Jo ^ 



(C+-e-^0(C_-e-^*; 

c+c_ 

C+C- 

(C+-e-^*)(C_-e-^*)' 



(A.2) 



(A.3) 



After executing the calculation, wc finally achieve 
(Nn) = {Nn)oe-'^^^^'''^^' 



(C+ -e-^*)(C_ -e-T*) 



+ 



(C+-e-^*)(C_-e-^* 



Lp[{l - p[)C+C_ 



1 _ 6-2(^^1) 



2(rT + ep;)-7 
2(rT + ep[) - 7 T 

27t _ -2{r^+ep{) 



+ ep[ - 7 



2L (1{C^ - p[{C. + C+)) - (C_ + C+)p;(l - Pi)) 
{C.-C+){Ng) 

(A.4) 



■?T + -(2p'i-l)+p'i(l-p'i))/> 
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for xi > X2, and 



(C+ -e-^*)(C_ -e-^*) 



+ 



1 _ e-2(ri-eP2) 



e-7« _ e-2(n-^P2) 27 



2(ri - ep'2) - 7 e 

e-27« _ g-2(r;-ep'2) 



+ 



(C_-C+)(7Vg) 



.2_ + 2(2p',_l)+p'^(l_p'^))L 



(A.5) 



for xi < X2. 
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